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Abstract 

We introduce the notion of globally updating evolution for a class of weighted networks, in 
which the weight of a link is characterized by the amount of data packet transport flowing through 
it. By noting that the packet transport over the network is determined nonlocally, this approach 
can explain the generic nonlinear scaling between the strength and the degree of a node. We 
demonstrate by a simple model that the strength-driven evolution scheme recently introduced can 
be generalized to a nonlinear preferential attachment rule, generating the power-law behaviors in 
degree and in strength simultaneously. 

PACS numbers: 05.40.-a, 87.23.Kg, 89.75.Hc 
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t research has arisen as the interdisciplinary subject for studying complex sys- 
DOGGO Although the binary (on/off) picture of connectivity has been 



shown to be quite informative and led to a number of important progress in our under- 
standing of complex systems, such as the ubiquity of power laws in its connectivity pattern, 
the degree distribution |8], Pd(k) ~ /c~ 7 , one may need to put a step forward to describe 
them more realistically. The weighted network, in which links between nodes in the net- 

'orm weights, is one of the most straightforward 
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may represent the intimacy between individuals in social networks, or the bandwidths of 
routers and optical cables in the Internet. The weighted network can be characterized by 
the generalized adjacency matrix W whose element {wij} denotes the weight of the link 
connecting the node i to j. By definition, Wij = when there is no link between i and j. 
We will restrict our interest here to the case of non-negative weight, {wsa > 0j for all (i, j) 
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(22J, |23J]. The strength s of a node given by 
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generalizes the concept of the degree, the number of links it has in binary networks. In 
terms of {?%}, the degree A; of a node may be written as 
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The analysis of weighted networks has been hindered mostly by the lack of large-scale 



data for the real networks. Recently, Barrat et al. [12| made the first detailed comparative 
analysis on the structure of weighted networks of the real world, the scientific coauthorship 
network (SCN) and the world-wide airport network (WAN). For the SCN, the weight of 
a link between two scientists is given roughly by the frequency of their collaboration, the 
effective number of papers they cowrote. For the WAN, the weight is taken as the total 
number of passengers of the direct flights between two connected cities. Interestingly, the 
two networks reveal qualitatively different organization of weight and network topology. For 
the SCN, the strength of a node (scientist) scales linearly with the degree, that is, 



s(k) ~ k . 



(3) 
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On the other hand, the WAN exhibits a nonlinear scaling as 

s(k) ~ k 13 , (4) 

with (3 Rj 1.5. 

Later, Barrat, Barthelemv. and Vespignani (BBV) introduced a simple evolntion model 
for such weighted networks |l4j . Basically, the BBV model is similar to the Barabasi-Albert 
(BA) model of binary scale-free (SF) network in spirit, containing the growth and the 
preferential attachment (PA) as basic ingredients. But the difference comes in the aspects 
that (i) strength at each node and weight at each link are introduced and (ii) the PA rule 
in the strength instead of the degree is applied. That is, the probability IT that an existing 
node i will receive a connection from a newly-introduced node is proportional to its strength, 
Ilj ~ Si. The strength of the target node subsequently increases by a constant amount A 
and so does the weight of the links incident upon the target node in a linear fashion. Then 
the degree and the strength scales linearly with each other and the distributions of them 
both follow power laws. On the other hand, the nonlinear scaling Eq. Q observed in the 
WAN requires different approaches. BBV proposed the nonlinear coupling between node- 
strength and link- weight, but a finite cutoff in node-strength was necessary to achieve the 
SF behavior in both the degree and the strength distribution j^J. 

In this Brief Report, we propose a packet transport driven evolution model of weighted 
SF network, which can illustrate the nonlinear relationship Eq. (4). We first point out that 



the weight used for the WAN in Ref. [12J is actually the amount of traffic through the 
link. From this perspective, the evolution of such weighted network should be viewed as 
packet transport driven. The distinguishing point here is that the traffic flowing through the 
network is determined in a nonlocal manner, in high contrast to the local weight evolution 
rule as formulated in the BBV mode l |l4l| . A measure of such traffic over the (binary) 



network is the quantity called the load [25], or the betweenness centrality |26|. The load of a 
node, the vertex-load, is defined by the effective number of data units passing through that 
node when every pair of nodes in the network sends and receives a unit data in unit time 
step. The data are assumed to be delivered only along the shortest path(s) between the 
pair. When the data encounters a branching point during the transport, they are supposed 
to be divided evenly by the number of branches. Thus the load quantifies the level of traffic, 
albeit in an idealized way. One can also define the link-load in a similar fashion, as the 
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FIG. 1: (a) The link-load vs. the product of the two degrees kikj of the nodes at each end of 
the link in the BA model with the system size N = 10 3 . Data points are logarithmically binned 
to reduce fluctuation. The straight line with slope 0.63 is drawn for reference, (b) The relation 
between the node strength (the sum of the link-loads attached to the node) and the load of the 
node. As indicated by the guideline with slope 1, they scale linearly for large I. 



effective number of data units passing through a given link. For SF networks, t 
each node is heterogeneous, and its distribution follows a power law, Pl(£) ~ C-~ 6 
the ranks of a node in degree and in load are preserved, then the scaling relation 
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holds for each node i, where rj = (7 — — 1). This relation is valid in the BA-type model 
and the Internet 2^ |. 

The weight of a link Wij and the product of the degrees of the nodes at its ends kikj are 
found to be related as a power law, Wij ~ (kikj) 9 , and 9 pa 1/2 for the WAN [l2j. Inter- 
estingly, the same half-power scaling was observed in the metabolic network of Escherichia 
coli, between the metabolic flux of a reaction and the degrees of the participating metabo- 
U to Q. It is — e to oote t h at Su ch a —oa has also been observed in the 
relation of the link- load versus kikj for binary networks 3> 29]. We show in Fig. 1(a) 
such a relation for the BA model Q], finding a slightly different scaling exponent 9 pa 0.63. 
So we may regard the link-load in the binary network as an approximation of the weight in 



the weighted networks characterized by the traffic level, such as the WAN. Although actual 
traffic level may well depend also on more complicated factors such as the geographic dis- 
tances between the airports and the queuing and transit delays at the airports, we exploit 
this idea as a starting point for further discussion below. In this setting, the strength of a 
node is given approximately by the load of the node itself, i.e., 

Si = y^Jij (6) 

since the vertex-load of a node is roughly one half of the sum of the link-loads of the links 
connected to that node for large i as shown in Fig. 1(b). 

When the relation (JSJ) holds, the strength and the degree of a given node scale nonlinearly 
as Si ~ k^. Then the application of BBV scheme leads to Hi ~ Sj ~ with 77 > 1 in most 
cases. Such a super-PA in degree breaks the power-law degree distribution |30J] as well as 
Eq. (0) itself. To obtain an SF network, one must use the sub-PA in strength, 

n, ~ sf ~ if (7) 

with a = I/77. At a first glance, one may think that this rule of sub-PA in strength would 
not generate a SF strength distribution, however, we can achieve a SF strength distribution 
with the strength being the load. 

To demonstrate the above idea, we simulate a growing network as follows, (i) Initially, 
mo = 3 nodes are introduced and they are fully connected, and we calculate the load at 
each node, (ii) At each time step, a new node is added, and attaches m = 2 links to existing 
nodes selected following the PA rule IT ~ if. (iii) The load of each node is recalculated. 
This process is repeated N times. Figs. 2(a-b) show the simulation results, the degree and 
the strength (or load) distributions of the packet transport-driven growth network for the 
cases of a,\ — 1 and a 2 = I/77 ~ 0.6. Indeed we can see that for the a = 1 case, the power- 
law behaviors of the degree and the strength distributions break down, while for a = 0.6, 
the model reproduces both the power laws with the degree exponent 7^3 and the load 
exponent 5 ~ 2.2, equivalent to those of the BA model. These measured values of 7 and 
5 are consistent with a = 0.6 through the relation a = (8 — l)/(7 — !)• Other structural 
properties such as the clustering and the degree mixing are also similar to those of the BA 
model [Figs. 2(c-d)]. It is of no surprise because this model with a = 1/rj is nothing but the 
BA model rephrased in terms of the strength-driven evolution. 
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FIG. 2: The cumulative distributions for the degree (a) and the strength (b) of the packet transport- 
driven network growth with Eq. (7). The open circles denote the data for a = 1 and the filled 
circles for a = 0.6. The guidelines have slopes —2.0 (a) and —1.2 (b). Note the humps for large 
degrees and large strengths in the case of a = 1, indicating the breakdown of the SF nature. 
The local clustering function C(k) (c) and the average neighbor degree function k nn (k) (d) for the 
network with a = 0.6. The filled circles denote the weighted version of the corresponding quantities 
introduced in Ref. and the cross symbols the original binary version of them. 



6 



Note that in our model, the link weight update is carried out through the change of the 
load, which occurs by global reorganization of shortest pathways within the network. As 
pointed out earlier, the load accounts for the traffic level only in an approximated manner. 
The general framework described in the paper, however, does not depend on what to use as 
the traffic measure as long as the scaling relation Eq.pjl holds for that quantity. 

Surely, it is not always the case that the weight of a link is determined by the level 
of traffic. In the SCN, for example, the weight measures the direct affinity between the 
scientists, which is primarily determined by their own attributes which is local in character. 
Indeed, for SCN we can see no appreciable correlation between the weight and the link-load 
of a link, as shown in Fig. |H1 In such a case, the BBV scheme could be more suitable than 
ours. Thus in the weighted network modeling, one has first to discriminate properly what 
the nature of the weight is in the system one wants to describe. 

Finally, we like to emphasize again the difference between the strength used here and that 
of the BBV model. We also compare them with the vertex-load directly measured. For the 
purpose, we measure the strength Sj of node i based on Eq. (1) where the weight Wij is given 
as the link- load £y, that is, Sj = an d compare it with the vertex-load £; measured 

directly from the BBV network and the quantity s^ BBV ^ defined in the BBV model [14] . We 
note that Si and ii (s\ ) are updated globally (locally) when a new node is added in the 
system. As shown in Fig. 01 we can find the scaling behaviors of s ~ i ~ A; 133 , which is 
consistent with the one obtained from the formula, I ~ ^(t- 1 )/^- 1 )^ anc i s ( BB v) ^ Those 
results support our claim. 

To conclude, we have introduced the notion of globally updating evolution in a class of 
weighted networks, in which the weight is characterized by the level of traffic flowing through 
the links. The link-load is used as a measure of the traffic in our model. This model explains 
the generic nonlinear scaling between the strength (traffic) and the degree, observed in, e.g., 
the WAN. We have also shown that the generalization of strength-driven evolution into the 
nonlinear PA in strength is necessary to produce a network which is SF both in degree and 
in strength. 

Note in added- After the first draft of the manuscript has been completed, we have learned 
of a recent work by Bianconi which introduced a model generating the nonlinear scaling 
relation between strength and degree, Eq.(£Q). The model succeeds to obtain the nonlinear 
relationship by decoupling the evolution of the network topology (degree) and the weights. 
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FIG. 3: Histogram of the links having the link load and the link weight in the coauthorship 
network of cond-mat subset of arXiv.org [3Jj. The grey level denotes the number of links in each 
bin, in logarithm with base 10. 

Although that model seems to be different from our conceptual model at first, the two 
models share important aspects. First its binary network structure grows following the BA 
model and second, the weights are updated nonlocally rather than locally as was assumed in 
the BBV model. Moreover, Guimera and Amaral [33[ introduced a model to illustrate the 
traffic in the WAN. The model also contains the rule of adding internal links between existing 
nodes, which are chosen by the PA rule combined with the ingredient of reducing traveling 
length. Such a rule is applied to all nodes in the network, not limited to neighbors of a 
newly added node, which needless to say accounts for global updating of traffic pathways. 
Such global reorganization of weights produces the non-linear relationship between strength 
and degree, which is the main conclusion of the work. 
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FIG. 4: The strength s(k) based on the link-load (□), the load i{k) based on the vertex-load (A), 
and the BBV strength s^ BY \k) (o) versus the degree k for the BBV model network with system 
size N = 10 4 , the number of links emanating from a new node m = 1, and the weight increment 
w o = A = 1.0. The slopes of the guide lines are 1.33 (dotted) and 1.0 (dashed), respectively. 
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